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AFRL  AURORAL  BOUNDARY  ALGORITHM 


1.  SUMMARY 

The  equatorward  auroral  boundary  location  is  determined  from  the  AFRL  algorithm  using 
electron  flux  data  obtained  with  the  SSJ5  sensor  onboard  DMSP  satellites.  This  report  presents 
the  data  used  for  determination,  describes  the  operation  of  the  algorithm,  and  describes  in  detail 
the  method  of  determination  of  the  four  equatorward  auroral  boundaries  for  orbit  9  of  DMSP 
satellite  FI 6,  25  July  2004  (Julian  day  207).  Results  are  presented  numerically  and  also  as 
spectral  views  of  these  boundary  locations  in  time  and  in  magnetic  latitude. 


2.  INTRODUCTION 

The  dynamics  of  the  auroral  oval  negatively  impact  the  operations  of  Air  Force  communications 
and  surveillance  systems  at  high  magnetic  latitudes.  Research  aimed  at  understanding  the 
vagaries  of  these  space  weather  effects  on  system  performance  led  to  the  development  of  new 
methods  to  specify  and  predict  the  location  and  intensity  of  auroral  optical-emission  and  plasma 
distributions.  By  the  early  1970s  auroral  scientists  recognized  that  the  nightside  auroral  oval 
magnetically  maps  to  the  magnetospheric  plasma  sheet.  Electrons  and  ions  with  energies  in  the 
kilo-electron  Volts  (keV)  range  in  the  plasma  sheet  precipitate  into  the  high-latitude  ionosphere 
to  create  optical  emissions  and  new  ionospheric  plasma. 

Using  auroral  electron  flux  measurements  by  the  particle  spectrometer  on  the  Defense 
Meteorological  Satellite  Program  (DMSP)  satellites  F2  and  F4,  Gussenhoven  et  al.  [1980,  1981] 
demonstrated  that  the  location  of  the  equatorward  boundary  of  the  auroral  oval  linearly  correlates 
with  the  Kp  geomagnetic  index.  The  prevailing  linear  relationships  vary  from  one  magnetic 
local  time  sector  to  the  next.  Gussenhoven  and  coworkers  at  AFRL  developed  automatic 
techniques  to  recognize  the  location  of  auroral  boundaries  from  variations  in  the  data  streams  of 
the  SSJ  sensors  that  fly  on  DMSP  satellites,  [ Hardy  et  al.,  1984].  The  algorithm’s  utility  is 
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limited  by  the  fact  that  boundary  identifications  are  observed  only  along  the  spacecraft’s 
trajectory.  Boundary  extensions  to  other  magnetic  local  times  derive  from  statistical  inference. 

DMSP  FI 6  was  launched  in  October  2003  carrying  the  Special  Sensor  Ultraviolet 
Spectrographic  Imager  (SSUSI).  In  principle  the  SSUSI  imager  offers  the  possibility  of 
extending  in  local  time  our  instantaneous  knowledge  of  the  auroral  boundary  location  via  direct 
UV  measurements.  To  validate  SSUSI’s  ability  to  identify  auroral  boundaries,  we  compared  its 
measurements  with  those  by  the  SSJ5  sensor  on  DMSP  FI  6.  In  support  of  this  effort  we  deemed 
it  useful  to  write  a  detailed  description  of  the  AFRL  algorithm  for  identifying  auroral  boundaries. 

The  AFRL  Auroral  Boundary  Algorithm  provides  automatic  identification  of  the  equatorward 
boundary  of  the  auroral  oval.  This  report  describes  (a)  the  database  used;  (b)  the  statistical 
parameters  generated  from  the  database;  (c)  the  actual  tests  performed  on  these  statistics;  (d)  the 
use  of  these  test  results.  Collectively,  these  details  constitute  the  operation  of  the  algorithm 
which  results  in  the  equatorward  auroral  boundary  determination.  The  specific  database 
provides  for  determination  of  4  boundaries  in  each  orbit  of  6060  seconds.  Boundaries  for  one 
complete  orbit  are  shown  as  they  are  generated,  together  with  the  database,  statistical  parameters, 
tests  results,  and  description  and  display  of  final  boundary  determination.  These  boundaries  are 
shown  also  graphically  as  spectral  images  in  both  compressed-time  and  expanded-time  views. 


3.  THE  DATABASE 

The  database  used  in  this  description  results  from  data  recorded  by  the  SSJ5  sensor  onboard  the 
DMSP  Satellite  FI 6,  during  orbit  9  of  the  satellite  on  25  July  2004  (Julian  Day  207).  The 
specific  data  used  for  boundary  determination  are  electron  counts  as  recorded  in  20  discrete 
energy  bins  of  the  J  sensor.  The  calculations  are  performed  by  considering  significant  changes  in 
slope  of  electron  count-rate  (counts/observation)  spectra,  a  process  dominated  by  data  acquired 
in  high-energy  channels. 
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4.  THE  SENSOR 

The  J  sensor  stores  electron  counts  in  20  assigned  energy  channels,  Ei  through  E2o,  where  the 
high-energy  channels  are  Ei  =  30keV,  E2  =  20keV,  E3  =  14.4keV,  E10  =  IkeV,  and  the  low- 

energy  channels  are  Eu  =  IkeV,...,  E2o=  30eV.  The  electron  counts,  ej (i  =  1 .. .20),  in  channel  Ei 
are  summed  every  4  seconds.  Times  assigned  to  the  sums  are  those  of  the  start  of  each  4-sec 
interval.  There  are  1515  such  4-sec  intervals  in  one  orbit  of  — 1 02  minutes,  and  14  orbits/day. 
The  counts  are  summed  in  pairs,  n,,  (i  =  1 ...  10),  where  ni  =  ei  +  e2,  n2  =  e3  +  e4,...,  nio  =  ei9+  e2o. 
For  DMSP  satellites  designated  FI 6  and  higher  there  is  no  channel  Eu-  Consequently,  in 
summing  the  counts,  since  n6  =  ei  1  +  ei2,  ei  1  has  been  set  equal  to  ej2. 

5.  THE  SATELLITE  ORBIT 

The  DMSP  satellite  orbit  is  considered  as  starting  at  the  equator  in  the  Evening  Sector  (see 
Figure  la).  As  the  satellite  moves  from  equator  to  pole  (poleward),  the  equatorward  edge  of  the 
auroral  boundary  is  approached.  Data  (electron  counts)  are  run,  as  collected,  through  the 
algorithm  in  quadrant  1  (Quad  1)  and  quadrant  3  (Quad  3)  (shaded  sectors).  Conversely,  in 
(Quad  2)  and  (Quad  4)  (white  sectors)  data  are  collected  as  the  satellite  moves  from  pole  to 
equator  (equatorward).  In  the  white  sectors,  the  sensor  encounters  auroral  activity  early  in  the 
sector,  and  exits  the  equatorward  boundary  later  in  the  sector  as  it  approaches  the  equator. 
Therefore,  for  such  sectors  the  data  are  run  through  the  algorithm  in  order  from  later-taken  to 
earlier-taken  (run  backwards  in  time  through  the  algorithm)  as  the  satellite  moves  equatorward. 
In  summary,  the  data  are  run  forward  in  the  shaded  sectors,  and  backwards  in  the  white  sectors. 
Each  of  the  4  boundary  calculations  uses  1515  seconds  of  data,  and  has  an  uncertainty  of  4 
seconds  (30  km)  for  J  sensors. 
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6.  STATISTICAL  PARAMETERS 

The  summed  counts,  ni,  become  the  basis  for  generating  a  set  of  five  statistical  parameters  Si, 
which  are  controlling  factors  in  the  boundary  determination  (see  Statistical  Parameter  Equations 
for  Si...  S5  below).  Each  S,-  is  tested  against  a  corresponding  threshold,  hj,  for  acceptance  or 
rejection.  The  thresholds  are  determined  for  the  J  sensor  on  each  satellite  during  post-launch 
calibration  in  the  laboratory.  In  this  report,  since  the  data  have  come  from  the  J5  instrument  on 
DMSP  satellite  FI 6,  the  five  thresholds  for  hi...  hs  are  350,  28,  20,  25,  and  35  respectively.  If  Si 
>  hi,  then  the  satellite  has  probably  encountered  a  boundary.  This  hypothesis  is  tested  to 
determine  whether  or  not  a  boundary  was  crossed  by  looking  for  consistency  in  exceeding  the 
threshold  value,  and  is  then  either  accepted  or  rejected  on  this  basis.  When  all  criteria  are 
satisfied  the  boundary  is  considered  crossed  and  the  boundary  search  is  completed. 

Note: 

The  summed  counts,  n„  and  the  statistical  parameters,  Si:  together  with  other  parameters  constitute  the  Database  for 
this  report,  and  are  shown  in  Table  2.  Due  to  printer  limitations  at  the  time  of  planning  and  testing  this  algorithm, 
where  Sj  number  sizes  become  too  large  to  fit  on  the  printed  page,  ****  appears  in  those  locations  in  this  Table. 
This  has  no  bearing  on  the  boundary  determination. 
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7.  STATISTICAL  PARAMETER  EQUATIONS 

Si  is  a  function  of  n^....  ng  (lower  energy  channel  counts),  but  the  other  four,  S2...  S5,  are  functions 
of  ni  ns  (higher  energy  channel  counts)  only.  The  following  equations  were  created  and 
chosen  strictly  by  trial  and  error  at  the  time  of  the  algorithm  development,  as  adjustments  were 
made  to  best  determine  the  boundary. 


Sj  =  n6  +  n7  +  «8  +  n9 


10  31m, -m,2  .  w,  =16/i,+8«2+4«3  +  2«4+/75 

—  - - - —  where  ,  ,  ,  ,  , 

31  J  ^1  YYi-,  —  1 6«,  +  8h2  +  +  2/?4  + 

V  31 


where 


m]  =  +  2«2  +  4«3  +  8h4  + 1 6n5 

m2  =  «,2  +  2  h22  +  4«32  +  8«42  + 16«52 


where 


w,  =  h,  +  n2 
m2=n4+  n5 
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8.  BOUNDARIES  AND  THEIR  DETERMINATION 


Boundaries  determined  from  the  data  recorded  by  DMSP  FI 6  on  25  July  2004  (Julian  day  207) 
are  shown  in  Table  1  for  the  entire  day.  The  bordered  area  encloses  four  boundaries  determined 
during  one  full  satellite  orbit.  The  values  for  these  boundaries  are  shown  also  in  Figure  lb 
below.  In  Quad  1  of  this  orbit  on  the  nightside,  termed  “evening  sector”  starting  at  47613  +  500 
sec  (48113  sec),  a  boundary  was  determined  at  48461  sec  at  55°. 3  Magnetic  Latitude,  157°. 9 
Magnetic  Longitude;  in  Quad  2  on  the  dayside,  termed  “morning  sector”  starting  at  49157  +  0 
sec  (49157  sec),  as  DMSP  moved  equatorward  from  the  pole  (N),  a  boundary  was  determined  at 
49693  sec  at  60°.  1  Magnetic  Latitude,  359°. 2  Magnetic  Longitude;  in  Quad  3,  on  the  dayside, 
“morning  sector”  starting  at  50673  +  600  sec  (51273  sec),  as  DMSP  moved  poleward  (S)  from 
the  equator,  a  boundary  was  determined  at  51617  sec  at  -56°. 0  Magnetic  Latitude,  338°. 6 
Magnetic  Longitude;  in  Quad  4,  on  the  nightside,  “evening  sector”  starting  at  52277  +  0  sec 
(52277  sec),  as  DMSP  moved  equatorward  from  the  pole  (S),  a  boundary  was  determined  at 
52833  sec  at  -57°. 0  Magnetic  Latitude,  150°. 6  Magnetic  Longitude. 


Quad  1 

Quad  2 

Quad  3 

47613  + 

!  49157  - 

50673-  1 

500 

Li  J 

600 

V 

Orbit  Start  »»> 

48113 

49157 

51273 

52277 

47613  + 

49157  + 

50673  + 

52277  4 

848 

536 

944 

556 

V 

V~ 

Boundary  - ► 

48461 

49693 

51617 

52833 

49157+ 

49157- 

52277- 

52277- 

0 

996 

0 

896 

V 

V 

Orbit  End  »»> 

49157 

50153 

52277 

53173 

Figure  lb 

Orbit  Sector  Start  Time,  End  Time,  and  Boundary  Time  as  Determined  for  Orbit  9 
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9.  BOUNDARY  ANALYSIS 

For  a  boundary  determination  to  occur,  there  must  be  significant  changes  in  the  slope  of  the 
count-rate  as  indicated  by  Sj.  The  data  used  for  the  determination  of  these  four  boundaries  are 
shown  on  the  Boundary  Analysis  sheets  of  Table  2  (21  pages),  including  n*  and  Sj  values  for 
every  4-sec  interval  of  the  orbit.  For  discussion  and  illustration  of  the  operation  of  the 
algorithm,  the  boundary  determination  in  Quad  4  (above)  will  be  examined  in  detail,  and 
revealed  by  shaded  and  colored  areas  on  page  31  (Table  2)  .  The  full  set  of  data  and  statistics 
for  this  boundary  determination  are  shown  on  pages  30-34  (Table  2). 

9.1  Boundary  Analysis  Table  Structure  -  Quarter  Orbits  Quad  1  and  Quad  3 

On  page  14,  the  first  page  of  the  Boundary  Analysis  Table  (Table  2),  the  25  data  columns  are 
labeled  ti...t25.  These  are  25  4-sec  intervals  of  time  for  listings  of  r\\  (i  =  1...9).  For  listing 
purposes,  recording  of  the  start  time  of  this  quarter  orbit  (Quad  1)  is  shown  as  47613  (sec).  On 
the  line  following,  observe  the  notation  “AT  T  =  500”  (sec).  Consequently,  the  time  for  entry  ni 
as  shown  at  ti  is  47613  +  500;  at  n2,  is  47613  +  504;  n3,  is  47613  +  508,  etc  across  the  row  of  25 
entries  -  time  increasing  by  4  sec  for  each  interval.  Therefore,  time  has  increased  by  100  sec  by 
the  start  of  the  data  on  the  lower  half  of  page  14  (“AT  T  =  600”).  This  proceeds  through  page  19 
of  Table  2,  to  the  end  of  Quad  1.  Recall,  from  Figure  la,  that  data  are  run  forward  through  the 
algorithm  in  Quad  1  and  Quad  3  (the  shaded  sectors).  Accordingly,  see  the  data  of  Quad  3 
(Table  2,  pages  25-29)  and  observe  that  time  increases  for  this  quarter  orbit  as  it  did  for  Quad  1 . 

9.2  Boundary  Analysis  Table  Structure  -  Quarter  Orbits  Quad  2  and  Quad  4 

Data  for  Quad  2  starts  on  page  20  of  the  Boundary  Analysis  Table  (Table  2).  For  listing 
purposes,  recording  of  the  start  time  of  this  quarter  orbit  is  shown  as  49157.  On  the  line 
following,  observe  the  notation  “AT  T  =  996”  (sec).  Consequently,  the  time  for  entry  ni  as 
shown  at  fi  is  49157  +  996;  at  n2,  49157  +  992;  at  n3,  49157  +  988,  etc  across  the  row  of  25 
entries  -  time  decreasing  by  4  sec  for  each  interval.  Therefore  time  has  decreased  by  100  sec  by 
the  start  of  the  data  on  the  lower  half  of  page  20  (“AT  T  =  896”).  Recall,  again  from  Figure  la, 
that  data  are  run  backwards  through  the  algorithm  in  Quad  2  and  Quad  4  (the  white  sectors). 
Accordingly,  view  the  data  of  Quad  4  (Table  2,  pages  30-34)  and  observe  that  time  decreases  for 
this  quarter  orbit  as  it  did  for  Quad  2. 
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10.  Quad  4  BOUNDARY  LOGIC  and  DETERMINATION 

Refer  to  page  30  (Table  2)  for  start  of  this  detailed  discussion  which  results  in  Quad  4  boundary 
determination  parameters  on  page  31.  Here  ti  =  52277  +  896  sec  [quarter  orbit  start  time  + 
number  of  seconds  elapsed  beyond  start  time  (Del  T)],  and  start  of  4-sec  summing  for  m;  t2  = 
52277  +  892  sec,  etc.  Special  note  must  be  made  here  of  page  31  where  the  boundary  will  be 
met,  and  the  lines  marked  Del  T,  Mlat,  S1...S5,  and  LOGIC.  Mlat  =  10  x  magnetic  latitude  in 
degrees  to  one  decimal  place.  This  form  of  Mlat  is  retained  in  Table  2  in  order  to  preserve  the 
format  of  some  printing  constraints  encountered  during  algorithm  development. 

10.1  LOGIC 

The  logic  state  (LOGIC)  for  any  given  4-sec  interval,  set  to  0  at  the  start  of  every  quarter  orbit,  is 
determined  by  the  Sj>  hjtest.  The  possible  LOGIC  for  a  4-sec  interval  is  0,  1,  -1,  1,  2,  or  3.  If  Sj 
>  hj,  q;  is  TRUE  for  that  Sj.  The  time  assigned  is  that  of  the  start  of  the  4-sec  interval.  If  Sj  >  hj 
for  3  consecutive  4-sec  intervals,  Q;  is  TRUE  for  that  Sj,  LOGIC  is  set  to  1,  and  the  time 
assigned  is  that  of  the  start  of  the  12-sec  interval.  That  LOGIC  never  returns  to  0  during  this 
quarter  orbit,  but  it  sets  to  -1  if  the  TRUE  condition  is  not  maintained  for  a  sufficient  duration. 

Possible  LOGIC  transitions 


0 

to  1 

One  or  two  QjS  are  TRUE 

1 

to  -1 

TRUE  for  QjS  does  not  maintain 

-1 

to  1 

One  or  two  Qi  s  are  TRUE  again 

1 

to  2 

Three  QjS  are  TRUE 

2 

to  3 

Five  Q;s  are  TRUE 

For  a  boundary  to  be  determined,  LOGIC  must  reach  2  or  3,  values  from  which  it  cannot  return 
to  a  lower  state.  When  LOGIC  =  2  is  reached,  three  QjS  are  TRUE.  While  in  LOGIC  =  2  state, 
if  data  are  exhausted  before  five  QjS  are  TRUE,  the  boundary  is  determined,  and  search  is 
complete.  To  reach  LOGIC  =  3,  Qi,  Q2,  Q3,  Q4,  and  Q5  must  all  be  TRUE,  the  boundary  is 
determined,  and  search  is  complete.  As  Qi  becomes  TRUE  for  its  Sj,  Magnetic  Latitude  (MLAT) 
of  the  boundary,  as  contributed  by  that  Si,  is  recorded  in  its  appropriate  position  (replacing  -999) 
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in  Table  2.  Flowcharts  of  this  operation,  leading  to  the  boundary  determination,  are  shown  in 
Figures  2  through  6  for  determination  of  Qi...  Qs  TRUE,  and  will  be  followed  for  here  for 
S1...S5  for  the  bottom  half  of  page  31  to  where  the  Evening  Sector  boundary  occurs,  for  shaded 
sector  Quad  4,  for  DMSP  FI 6  Orbit  9,  Julian  Day  207  (25  July  2004). 

10.2  DETERMINATION 
Observations 

In  the  top  half  of  page  30,  within  the  black-bordered  region  (where  Sj  values  are  enclosed),  the  Si 
>  hj  test  has  never  passed  for  a  4-sec  interval,  the  qi  has  never  been  TRUE,  and  certainly  S;  >  h; 
for  three  consecutive  4-sec  intervals  has  not  been  TRUE.  Therefore,  Qi  has  not  been  TRUE  for 
any  Sj.  The  same  situation  holds  for  the  bottom  half  of  page  30.  Everything  of  interest  to  us  in 
this  discussion  happens  on  page  31 .  Refer  to  Figures  2  -  6  for  flowcharts  generating  +,  *,  and  A 
as  discussed  below. 

Si  >  Test  Results  and  Analysis 

The  first  incidence  of  Si  >  h,  TRUE  in  this  quarter  orbit  occurs  for  S2  at  Del  T  =  656,  and  is 
reflected  by  a  “+”  after  the  S2  value  of  30.  At  the  same  Del  T,  Si  >  hi  TRUE  for  S3,  as  evidenced 
by  a  “+”  after  the  S3  value  of  25.  In  the  case  of  S2,  the  Sj  >  hj  TRUE  does  not  maintain  at  Del  T  = 
652,  the  next  4-sec  interval  (as  seen  by  a  lack  of  after  the  S2  value  of  25).  Observe  the 
behavior  of  S2,  S3,  and  S4  for  the  remainder  of  the  top  half  of  page  31.  Q1...Q5  have  never 
proven  TRUE  for  the  span  of  time  from  Del  T  =  696  through  Del  T  =  600.  Note  that  for  this 
case  if  S2  >  I12  TRUE  had  maintained  for  three  consecutive  4-second  intervals,  the  three  cells 
would  have  contained,  in  sequence,  a  “+”  a  a  “A”  after  the  S2  value  in  each  respective  cell. 

On  the  bottom  half  of  page  31,  activity  increases  greatly.  Note  that  LOGIC  has  gone  from  0  to  -1 
in  column  2  to  column  3.  This,  however,  is  not  an  allowed  transition.  Transition  from  0  to  1  is 
allowed,  never  returning  to  0  again  in  any  given  quarter  orbit.  Also,  transition  1  to  -1  is  allowed 
when  the  TRUE  condition  does  not  maintain.  What  has  happened  here  is  that  transition  0  to  1  to 
-1  has  occurred  too  rapidly  to  register  the  0  to  1  transition.  Since  S2...S5  are  have  not  changed 
considerably  in  this  time  frame,  the  0  to  -1  transition  most  probably  had  occurred  due  to  change 
in  the  Si  value. 


9 


Now  note  the  following: 


at  Del  T  =  560, 

at  Del  T  =  556, 

at  Del  T  =  552, 

at  Del  T  =  552, 

at  Del  T  =  528, 


54  starts  a  sequence  of  “+”  “A” 

51  starts  a  sequence  of  “+”  “A” 

52  starts  a  sequence  of  “+”  “A” 

55  starts  a  sequence  of  “+”  “A” 

53  starts  a  sequence  of  “+”  “A” 


The  first  full  sequence  is  completed  by  S4  at  Del  T  =  552,  and  so  LOGIC  transitioned  from  -1  to 
1  at  that  time.  At  that  time  Si,  S2,  and  S5  began  or  continued  full  sequences,  and  LOGIC 
remained  at  1.  When,  finally,  S3  completed  a  sequence  at  Del  T  =  520,  LOGIC  transitioned  to  2 
at  that  time,  since  now  Qi,  Q2,  and  Q3  are  TRUE.  Recall  from  the  Boundary  Analysis  section, 
that  although  the  LOGIC  is  set  at  the  times  reflected  on  page  31,  the  time  assigned  to  these  Qi 
TRUE  events  are  the  times  of  the  start  of  the  completed  12-sec  intervals. 


Boundary  Time  and  Magnetic  Latitude 

For  this  boundary,  the  Magnetic  Latitude  (Mlat)  for  times  for  Qi  TRUE  are  (yellow): 


Qi 

q2 

q3 

Q4 

Qs 

Mlat »» 

570 

572 

585 

568 

572 

Replacing»> 

-999 

-999 

-999 

-999 

-999 

(in  the  appropriate  row  on  page  31) 
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Steps  in  calculation: 

1 )  Average  the  Mlat  values  (MlatAvg) 

2)  Take  differences  (Mlat  -  MlatAvg) 

3)  If  all  within  1 0  difference,  assign  MlatAvg  as  the  Boundary  Magnetic  Latitude 

4)  If  not  all  within  1°  difference,  discard  Mlat  farthest  from  MlatAvgand  repeat  from  1) 

until  3)  is  satisfied. 


Calculation  Illustration 


Differences 

Differences 

Original 

Remaining 

Mlats 

5  entries 

Mlats 

4  entries 

56.8 

-0.5 

56.8 

-0.5 

57.0 

-0.3 

57.0 

-0.3 

57.2 

-0.1 

57.2 

-0.1 

57.2 

-0.1 

57.2 

-0.1 

58.5 

1.2 

58.5  eliminated 

286.7 

228.2 

Mlat  Avg»> 

57.3 

57.1 

#  Entries  »» 

5.0 

4.0 

Resulting  Mlat  is  truncated  to  57.0 


Time  corresponding  to  the  Boundary  Mlat  is  looked  up  in  the  ephemeris. 

Here,  page  31,  the  Del  T  is  seen  to  be  556.  The  Boundary  is  at  52883  seconds  (52277  + 
556). 


11.  SPECTRAL  IMAGES 

Spectral  images  of  the  four  orbit  9  boundary  determinations  (upper  portion)  and 
expanded  view  of  each  boundary  (lower  portion)  are  shown  in  Figures  7,  7a,  and  7b  as 
functions  of  UT,  Local  time,  and  Magnetic  Latitude. 


Figure  7  Compressed-time  view  of  all  4  boundaries 

Figure  7a  Expanded-time  view  -  Quad  1  and  Quad  2  boundaries 

Figure  7b  Expanded-time  view  -  Quad  3  and  Quad  4  boundaries 
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Table  1.  Auroral  Boundaries  for  entire  day  207 


Sat  ID 

Index 

UT  sec 

UT  hms 

Year 

Day 

Glat 

Glon 

Ml  at 

Mlon 

FI  6 

l 

449 

0:07:29 

2004 

207 

65.1 

139.0 

60.4 

206.8 

F16 

7 

2293 

0:38:13 

2004 

207 

-41.8 

103.7 

-57.0 

169.3 

FI  6 

9 

3453 

0:57:33 

2004 

207 

-68.0 

309.6 

-56.1 

16.6 

F16 

16 

5441 

1:30:41 

2004 

207 

46.6 

269.0 

59.4 

337.7 

F16 

17 

6557 

1:49:17 

2004 

207 

65.5 

114.0 

61.2 

185.3 

F16 

23 

8385 

2:19:45 

2004 

207 

-40.4 

78.7 

-53.7 

139.2 

F16 

25 

9613 

2:40:13 

2004 

207 

-65.7 

281.3 

-53.8 

0.4 

F16 

30 

11581 

3:13:01 

2004 

207 

48.0 

243.1 

57.0 

303.5 

F16 

31 

12697 

3:31:37 

2004 

207 

64.2 

87.2 

60.6 

160.9 

F16 

37 

14601 

4:03:21 

2004 

207 

-46.1 

51.1 

-55.0 

106.8 

F16 

40 

15837 

4:23:57 

2004 

207 

-59.8 

250.6 

-51.9 

338.3 

F16 

44 

17777 

4:56:17 

2004 

207 

52.5 

215.4 

55.7 

272.9 

F16 

46 

18853 

5:14:13 

2004 

207 

62.0 

59.8 

59.2 

135.1 

F16 

53 

20813 

5:46:53 

2004 

207 

-51.6 

23  2 

-55.2 

76.9 

F16 

55 

22077 

6:07:57 

2004 

207 

-52.8 

221.0 

-51.4 

310.0 

FI  6 

58 

24017 

6:40:17 

2004 

207 

59.5 

185.8 

57.5 

245.0 

F16 

59 

25005 

6:56:45 

2004 

207 

60.0 

32.7 

58.0 

110.4 

F16 

66 

27005 

7:30:05 

2004 

207 

-55.9 

355.4 

-53.7 

50.1 

F16 

68 

28269 

7:51:09 

2004 

207 

-48.5 

193.6 

-53.2 

281.3 

F16 

73 

30197 

8:23:17 

2004 

207 

63.0 

157.5 

58.8 

222.1 

F16 

74 

31197 

8:39:57 

2004 

207 

55.8 

4.5 

55.2 

84.9 

F16 

81 

33177 

9:12:57 

2004 

207 

-59.0 

328.0 

-51.0 

28.2 

F16 

83 

34477 

9:34:37 

2004 

207 

-43.2 

166.0 

-53.7 

248.7 

F16 

87 

36257 

10:04:17 

2004 

207 

59.9 

134.5 

55.6 

204.2 

F16 

89 

37325 

10:22:05 

2004 

207 

55.1 

338.6 

57.7 

63.4 

F16 

96 

39381 

10:56:21 

2004 

207 

-63.9 

298.6 

-52.1 

9.9 

F16 

98 

40649 

11:17:29 

2004 

207 

-40.0 

139.3 

-54.2 

215.5 

F16 

102 

42325 

11:45:25 

2004 

207 

57.2 

110.8 

53.7 

183.1 

F16 

104 

43437 

12:03:57 

2004 

207 

55.4 

313.3 

62.8 

39.2 

F16 

111 

45593 

12:39:53 

2004 

207 

-69.0 

267.3 

-57.8 

353.7 

F16 

113 

46813 

13:00:13 

2004 

207 

-37.3 

112.9 

-52.7 

182.2 

F16 

118 

48461 

13:27:41 

2004 

207 

58.3 

84.7 

55.3 

157.9 

F16 

120 

49693 

13:48:13 

2004 

207 

47.4 

283.9 

60.1 

359.2 

F16 

127 

51617 

14:20:17 

2004 

207 

-64.1 

247.5 

-56.0 

338.6 

F16 

129 

52833 

14:40:33 

2004 

207 

-42.8 

89.4 

-57.0 

150.6 

FI  6 

135 

54605 

15:10:05 

2004 

207 

59.9 

58.1 

57.3 

133.1 

F16 

136 

55885 

15:31:25 

2004 

207 

43.0 

256.7 

54.7 

321.5 

F16 

143 

57657 

16:00:57 

2004 

207 

-59.8 

225.5 

-56.4 

319.1 

F16 

145 

58913 

16:21:53 

2004 

207 

-44.9 

64.6 

-55.8 

121.2 

F16 

151 

60673 

16:51:13 

2004 

207 

57.1 

34.4 

55.2 

110.8 

F16 

153 

61937 

17:12:17 

2004 

207 

46.7 

232.7 

53.8 

292.4 

F16 

159 

63705 

17:41:45 

2004 

207 

-56.0 

202.5 

-58.0 

295.8 

F16 

161 

64945 

18:02:25 

2004 

207 

-49.7 

41.2 

-56.3 

94.8 

F16 

166 

66773 

18:32:53 

2004 

207 

56.2 

9.5 

55.3 

89.2 

F16 

168 

67913 

18:51:53 

2004 

207 

54.7 

211.0 

56.9 

267.9 

F16 

174 

69749 

19:22:29 

2004 

207 

-51.9 

179.2 

-59.3 

269.0 

FI  6 

175 

70929 

19:42:09 

2004 

207 

-57.2 

19.6 

-58.1 

69.9 

F16 

182 

72877 

20:14:37 

2004 

207 

55.5 

344.4 

57.2 

68.4 

F16 

184 

73957 

20:32:37 

2004 

207 

58.8 

188.0 

57.1 

247.0 

F16 

189 

75689 

21:01:29 

2004 

207 

-41.8 

157.8 

-53.7 

238.7 

F16 

191 

76949 

21:22:29 

2004 

207 

-62.5 

358.0 

-58.0 

49.2 

F16 

198 

78989 

21:56:29 

2004 

207 

55.3 

319.1 

61.7 

45.3 

F16 

200 

80001 

22:13:21 

2004 

207 

62.9 

165.6 

58.9 

228.3 

F16 

205 

81793 

22:43:13 

2004 

207 

-41.2 

132.6 

-55.9 

207.1 

F16 

207 

83021 

23:03:41 

2004 

207 

-65.0 

334.8 

-56.2 

32.0 

F16 

214 

84953 

23:35:53 

2004 

207 

46.6 

297.7 

57.7 

18.2 

F16 

216 

86085 

23:54:45 

2004 

207 

64.6 

141.8 

60.0 

209.2 
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Figure  2 

Flow  Chart  for  Sj  Boundary  Test 
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Figure  3 

Flow  Chart  for  S2  Boundary  Test 
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Figure  4 

Flow  Chart  for  S3  Boundary  Test 
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c  =  0  (counter  for  3  consecutive  Si  >  hi  values) 
Replaced  _ 
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Figure  5 

Flow  Chart  for  S4  Boundary  Test 
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Figure  6 

Flow  Chart  for  S5  Boundary  Test 
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Figure  7 

Compressed  Time  View  -  All  4  Boundaries 
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Figure  7a 


Expanded  Time  View  -  Quad  1  and  Quad  2  Boundaries 
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Figure  7b 

Expanded  Time  View  -  Quad  3  and  Quad  4  Boundaries 


40 


REFERENCES 


Gussenhoven,  M.  S.,  D.  A.  Hardy,  and  W.  J.  Burke,  Comment  on  "Diurnal  variation  of  the 
auroral  oval  size",  by  C.  -I.  Meng,  J.  Geophys.  Res.,  85,  2373,  1980. 

Gussenhoven,  M.  S.,  D.  A.  Hardy,  and  W.  J.  Burke,  DMSP  F2  electron  observations  of 
equatorward  auroral  boundaries  and  their  relationship  to  magnetospheric  electric  fields,  J. 
Geophys.  Res.,  86,  768, 1981. 

Hardy,  D.  A.,  L.  K.  Schmidt,  M.  S.  Gussenhoven,  F.  J.  Marshall,  H.  C.  Yeh,  T.  L.  Shumaker,  A. 
Huber,  and  J.  Pantazis,  Precipitating  electron  and  ion  detectors  (SSJ/4)  for  the  block  5 D/Flights 
4-10  DMSP  satellites:  Calibration  and  data  presentation.  Rep.  AFGL-TR-84-03 17,  Air  Force 
Geophys.  Lab.,  Hanscom  Air  Force  Base,  Mass.,  1984. 


41 


